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ABSTRACT 

We have determined the angular diameters of two metal-poor stars, HD 122563 and Gmb 1830, using CHARA and Palomar Testbed 
Interferometer observations. For the giant star HD 122563, we derive an angular diameter #30 = 0.940 ± 0.01 1 milliarcseconds (mas) 
using limb-darkening from 3D convection simulations and for the dwarf star Gmb 1830 (HD 103095) we obtain a ID limb-darkened 
angular diameter 6m = 0.679 ± 0.007 mas. Coupling the angular diameters with photometry yields effective temperatures with 
precisions better than 55 K (r eff = 4598 ± 41 K and 4818 ± 54 K — for the giant and the dwarf star, respectively). Including their 
distances results in very well-determined luminosities and radii (L = 230 ± 6 L Q , R = 23.9 ± 1.9 Rq and L = 0.213 ± 0.002 L Q , 
R = 0.664 ± 0.015 Rq, respectively). We used the CESAM2k stellar structure and evolution code in order to produce models that 
fit the observational data. We found values of the mixing-length parameter a (which describes ID convection) that depend on the 
mass of the star. The masses were determined from the models with precisions of <3% and with the well-measured radii excellent 
constraints on the surface gravity are obtained (logg = 1.60 ± 0.04,4.59 ± 0.02 dex, respectively). The very small errors on both 
logg and r e (f provide stringent constraints for spectroscopic analyses given the sensitivity of abundances to both of these values. The 
precise determination of r cff for the two stars brings into question the photometric scales for metal-poor stars. 

Key words. Stars: fundamental properties — Stars: individual HD 122563, HD 103095 (Gmb 1830) — Stars: low-mass — Stars: 
Population II — Galaxy: halo — Techniques: interferometry 



1. Introduction 

Metal-poor stars are some of the oldest stars in the Galaxy and 
thus reflect the chemical composition of Galactic matter at the 
early stages of Galactic evolution. The determination of accu- 
rate observed fundamental properties, and in particular their lo- 
cation in the Hertzsprung-Russell (HR) diagram, is a key re- 
quirement if we aim to constrain the unobservable properties 
such as mass, age, and initial helium content by using stellar 
models. Among the most controversial observed parameter is 
the effective temperature (T e s) which can vary by more than 
200 K for metal-poor stars from one meth od to another (see 
the PASTEL catalogue. lSoubiran et alj2010l) . In particular, local 
thermodynamic equilibrium (LTE) is usually assumed and non- 
LTE (NLTE) effects must be included in spectroscopic analyses 
especially for meta l -poor s tars where these effects a re enhanced 
(lThevenin&Idiarllll999r. lAndrievskv eTal] l2010b iMerle et alJ 



l20Tl and this leads to even more discrepancy between litera- 
ture values. One solution is to measure the angular diameter and 
convert this to T e ff to provide a direct determination. 

The large majority of metal-poor stars belong to the halo or 
the old disk of the Galaxy which means that their apparent mag- 
nitude and or angular diameters are extremely small and difficult 
to measure. Howev er, some instruments, in par ticular those on 
the CHARA array dten Brummelaar et al. 1 I2005I) are very capa- 
ble of working at short wavelengths on long baselines to ob- 
tain the required angular resolution. Among the most exciting 
possible targets with CHARA working in the K band are HD 
122563 (= HR 5270, HIP 68594, m v = 6. 19 mag) and Gmb 1830 
(= HD 103095, LHS 44, HIP 57939, m v = 6.45 mag) whose 
mean metallicities [Z/X]£] are — 2.3 dex and -1.3 dex, respec- 
tively (see discussion in Sect. 14. Il l, where Z and X denote the 



1 [Z/X] = logZ/X star - logZ/X and Z/X G = 0.0245, see Sect.l4~2l 
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Table 1. Most recent photometric and spectroscopic determina- 
tions of atmospheric parameters for the target stars 



Table 2. CHARA observations of HD 122563 and Gmb 1830. 
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^NLTE analysis by [M ashonki na et al.l J2008h Cuck & Heiteil d2006 
iMishenina & Kovtvukhl d2001h iGehren et all J20061) 



metallicity and hydrogen (absolute) mass fraction in the star and 
the subscript refers to the observed surface value. 

H D 122563, a standard example of a very m etal-poor field 
giant ( IWallerstein et al.l 119631: IWolfframl |1972|) . has been ex- 
tensively studied and presents similarities with metal-poor gi- 
ants found in globular clusters. Gmb 1830 is a metal-poor halo 
dwarf star recogni zed as exhibiting depleted Li dDelivannis et al.l 
Il994t I Kind 1 19971) when compared to the mea n value of halo 
dwarf stars dSpite & Spite] 119931: iRvanl 120051) . It is also the 
nearest halo dwarf and has an excellent parallax measurement. 
Combining interferometric measurements of these stars with 
other already measured old moderately metal-po or stars, such 
as // Cas ([Z/X] s = -0.5 dex. lBovaiian et aill2008l) . offers an ex- 
cellent opportunity to constrain the T e ff scale of metal-poor stars 
over a wide range of metallicities with possible implications for 
T e ff calibrations of globular cluster stars. In Table [TJ we summa- 
rize some of the most recent determinations of the atmospheric 
properties of both targets. Note that HD 122563 and Gmb 1830 
have also been defined as benchmark stars for the Gaia mission 
under the working group. 

Not only are temperature scales for metal-poor stars contro- 
versial, but stellar structure and evolution models often predict 
higher T e g than those observed for these stars (see e.g. Fig. 2 of 
lLebretorj2 000). The difficulty encountered when trying to match 
evolutionary tracks to the observational data not only severely 
inhibits the determination of any fundamental properties but any 
chance of improving or testing the physics in the models is also 
limited. 

Considering the difficulties mentioned above, in this pa- 
per we aim to determine accurate fundamental properties of 
HD 122563 and Gmb 1830 based on interferometric observa- 
tions (Sect. [2j. In Sect.[3]we present our analysis of the observa- 
tions to determine the angular diameters of both stars. We then 
determine the observed values of T e g, luminosity L, and radius 
R, and subsequently use stellar models to constrain the unob- 
servable properties of mass M, initial metal and helium content 
Z p Yi, mixing-length parameter a and age (Sect.|4|i. We also pre- 
dict their global asteroseismic properties in order to determine if 
such observations could further constrain the models. 



2. Observations 

The observations we r e col lected at the CHARA Array 
dten Brummelaar et"aT1 l200l . located at Mount Wilson 
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www. anst .uu. se/ulhei458/GaiaSAM/ 



Notes. MJD is the average modified julian date of the exposures and 
Inst, the instrument code (F: FLUOR, C: Classic). 



Observatory (California), together with two beam combining 
instruments: CHARA Classi c and FLUOR. CHARA Classic 
dten Brummelaar et al.1 [2005) is a two-telescope, pupil-plane, 
open-air beam combiner working in both the H and K' bands, 
and our observations correspond to the K' band (the centra l 
wavelength is A — 2.141 jum, from iBowsher et aD 1201 Oh . 
The raw data were reduced using th e pipeline described in 
ten B rummelaar et al. d200l . FLUOR dCoude du Foresto et al.l 
1 19981: iMerand et al.l I2006IT is a two-telescope beam combiner, 
but uses single-mode optical fibers for recombination. Single- 
mode fibers efficiently reduce the perturbations induced by 
the turbulent atmosphere on the stellar light wavefront, as 
the inject ed light corresponds only to the mo de guided by 
the fiber dRuilieri \l99^ ICoude du Forestol 119981) . Most of the 
atmospherically corrupted part of the wavefront is lost into the 
cladding, and the beam combination therefore occurs between 
two almost coherent beams. This results in an improved stability 
of the m easured fringe contr ast. The FLUOR data reduction 
pipeline (Me rand et al.l I2006L see also iKervella et all l2004bl) 
is based on the Fourier alg orithm and was developed by 
ICoude du Foresto et al.l d 19971) . 
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Table 3. FLUOR calibrator stars. 
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Table 4. CHARA Classic calibrator Stars 
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Sp. type 
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HD 119550 


G2V 


6.9 


5.3 


0.389 ± 0.027 


HD 122563 


HD 120066 


GOV 


6.3 


4.9 


0.479 ± 0.033 


HD 122563 


HD 120934 


A1V 


6.1 


6.0 


0.198 ±0.014 


HD 122563 


HD 121560 


F6V 


6.2 


4.8 


0.460 ± 0.030 


HD 122563 


HD 122365 


A2V 


6.0 


5.7 


0.238 ±0.016 


HD 122563 


HD 103799 


F6V 


6.6 


5.3 


0.343 ±0.013 


Gmb 1830 



We observed HD 122563 and Gmb 1830 in late 2007 and 
2008 with FLUOR and Classic. The corresponding visibility 
measurements V and uncertainties cr(V) are listed in Table [2] 
along with the projected baseline B and the baseline position 
angle PA measured clockwise from North. To monitor the inter- 
ferometric transfer function, we interspersed the observations of 
our two science targets with calibrator stars. The calibrators for 
the FLUOR observa tions were selected from the catalogue by 
iMerand et alj d2005l) . and these are listed in Table [3] and those 
for the CHARA Classic observations used the calibrators pre- 
sented in Table |U 

We also retrieved archival observations of HD 122563 in the 
K band obtained with the Palomar Testbed Interferometer (PTI) 
dColavitaetalJll999h between 1999 and 2002, and these are 
listed in Table [5] The data processing algorithm that was em- 
ployed Uore^luce thePTI observations has been described in de- 
tail bv lColavital (Il999l) . Due to the shorter baselines, the PTI ob- 
servations resolve HD 122563 marginally, and therefore do not 
strongly constrain its angular diameter. However, thanks to the 
relatively large number of observations, they provide an inde- 
pendent method for testing any bias in the CHARA observations. 



3. From visibilities to limb-darkened angular 
diameters 

We emp loyed a non-linea r, least-squares fitting routine in IDL 
(MPFIT. lMarkwardtl2009l) to fit uniform disk and limb-darkened 
visi bility functions for a single s t ar to the calibrated d ata points 
(see lHanburv Brown et alJll974l iBovaiian et al ] l2012h . We ob- 
tained a uniform disk diameter for HD 122563 and Gmb 1830 
of UD = 0.924 ± 0.01 1 mas and 6> UD = 0.664 + 0.015 mas, re- 
spectively. W e used the linear limb-darkening coefficients from 
lOarel d2000h assuming a [Fe/H]= -2.5, T eft - = 4500 K, and 
logg = 1.0 for HD 122563 and [Fe/H]= -1.5, T eS = 5000 K, 
and \ogg = 4.5 for Gmb 1830. The assumptions on these pa- 
rameters on the adopted coefficients have minimal influence on 
the final limb-darkened diameter, adding uncertainties of only 
a few tenths of a percent, well within the errors of our di- 
ameter measurements. We obtained #ld = 0.948 + 0.012 and 
6> LD = 0.679 + 0.0 15 for HD 122563 and Gmb 1830, respectively. 

We obtained a reduced^ 2 value of 0.28 for HD 122563 and 
0. 18 for Gmb 1830 from the fits. These values, much less than 1, 
are indicative of our individual measurement errors being over 
estimated. We show the data and the visibility function fits for 
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Fig. 1. Calibrated observations for PTI (black crosses), CHARA 
Classic (blue circles) and CHARA FLUOR (red squares) data 
plotted with the ID limb-darkened visibility function fit for 
HD 122563. See Section[3]for details. 
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Fig. 2. Calibrated observations for CHARA Classic (blue cir- 
cles) and CHARA FLUOR (red squares) data plotted with 
the limb-darkened visibility function fit for Gmb 1830. See 
Section[3]for details. 



HD 122563 and Gmb 1830 in Figures Q] and [2] The results from 
the fits to the data yield ID limb-darkened angular diameters of 
HD122563 and Gmb 1830 with a precision of -2% (see Table©, 
respectively. 

3.1. 3D limb-darkened angular diameter for HD 122563 

Convection plays a very important role in the determination of 
stellar limb-darkening. It has been shown that a 3D hydrody- 
namical treatment of the surface layers can lead to a significa- 
tive change of temperature gradients compared to ID hydrostatic 



sitv variation (e.g. lAllende Prieto et a 


. 20021 


Bieot et alj 


2006; 


Pereira et al. 


20091 IChiavassa et alj 


:oi(t 


Biaot et al.l 


2011; 


Havek et al.l 


2012). The 3D/1D limb-darkening correction for a 



I2OT0I) and is generally much stronger than for a dwarf star. We 
therefore used a radiative-hydrodynamical (RHD) surface con- 



3 



Creevey et al,; Stellar properties of metal-poor stars from interferometry 



Table 5. PTI observations of HD 122563. 
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vection simulation of a red giant for HD 122563 to determine 
the 3D limb-darkened angular diameter. The parameters of the 
model are (T e g) = 4627 + 14K (temporal average and standard 
deviati on of the effective temperature), |"Fe/H"|= -3.0, and logg 
= 1.6 dCollet et al.ll2009t IChiavassa et al.ll20Tob . The computa- 
tional domain of the RHD simulation represents only a small 
portion of the stellar surface (~l/30 of the circumspherence), 
however, it is sufficiently large to contain enough granules (~10- 
15) at each time step of the simulation. Hydrodynamical equa- 
tions are solved on a staggered mesh with a conservative schem e. 
Details of the computation can be found in lCollet et al.l (120091) . 

We computed emergent intensity for a representative se- 
ries of simulated snapshots and for wavelengths correspond- 
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Fig. 3. Best matching 3D-RHD synthetic visibility curves and 
PTI (black crosses), Classic (blue circles), and FLUOR (red 
squares) data for HD 122563. 



ing to the FLUOR filter (2.14 + 0.26 urn, equivalent to 
that for CHA RA) using the 3D pur e-LTBj radiative transfer 
code Optim3D (IChiavassa et alj|2009l) . It considers the Doppler 
shifts due to convective motions. Radiative transfer is solved 
monochromatically using pre-tabulated extinction coefficients 
for the same chemical compositions as the RHD simulations. It 
also uses the same extensive atomic and mole cular opacity data 
as the latest generation of MARCS models dGustafsson et al.l 
l2008h . 

For each time-step, we solve the radiative transfer equa- 
tion for different inclinations with respect to the vertical whose 
cosines are 1.000, 0.989, 0.978, 0.946, 0.913, 0.861, 0.809, 
0.739, 0.669, 0.584, 0.500, 0.404, 0.309, 0.206, 0.104]. From 
these limb-darkened intensities, we derived the monochromatic 
visibility curves using the Hankel Transform. The visibilities are 
then averaged with the transmission function of the instrument in 
the considered filter waveleng th domain. The pro cedure used in 
this work is the same as that of lBigot et al.l (1201 ll) . The synthetic 
visibilities are used to fit the interferometric K band observations 
given in Tables [2] and [5] 

Figure [3] displays the best fit of the visibility curve to the 
data that results in an angular diameter of #3d=0.940 ± 0.011 
mas (Tabled with a^ 2 = 0.35. Its value lies between that of the 
uniformed disk and ID limb-darkened diameters. This is a con- 
sequence of the fact that in realistic 3D hydrodynamical treat- 
ment of the stellar surface, the emergent intensity is less limb- 
darkened than the ID hydrostatic case. We note that the choice 
of the exact fundamental parameters of the 3D simulation does 
not influence the limb-darkened intensity and the derived angu- 
lar diameter by much. 

The 3D/1D correction is important f or determining the zero 
point of the effective temperature scale: IChiavassa et al.l (1201 Ol) 
(Table 3) showed that, in the case of metal-poor stars like the one 
analyzed in this work, 03d/#id ~ 2% in the K band (3.5% in the 
visible). This can result in corrections to the effective tempera- 
ture of ~40 K in the K band. In this case the resulting correction 
to the effective te mperature is ~15 K (see Sect.l4.11>. 

We note that iGonzalez Hernandez & Bonifa cio (2009]) pre- 
dict ID limb-darkened angular diameters for both of these stars 



3 Pure-LTE refers to when the source function is equal to the Planck 
function 
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Table 6. Angular diameters 





#of 


0VD-- 


t <T 


6iD-- 


t <T 


d 3D ±cr 


Star 


Observations 


(mas) 


(mas) 


(mas) 


Gmbl830 


18 


0.664 ± 


0.015 


0.679 ± 


0.015 




HD 122563 


66 


0.924 ± 


0.011 


0.948 ± 


0.012 


0.940 ±0.011 



using the infra-red flux method (IRFMfl For HD 122563 they 
predict 6»i D = 0.84 + 0.04 mas and for Gmbl830 lD = 0.61 
± 0.02 mas; both values are lower than our derived values. For 
HD 122563 this could be due to the fact that they use 2MASS 
pho tometry that is sa turated for this star (saturation limit K s ~ 

4.0. 1CutrietaIll200l . 

4. Constraints on stellar evolutionary models 

4.1. Observed parameters 

The list of observed parameters are summarized in the top 
part of Table 7] The magnitudes in the V band are take n from 
iJohnson et al. (Il966l) . those in the K b ands are from iDucatil 
(I2002I) for HP 122563 and ICutri et al J (120031) fo r Gmb 1830, 
and the Hipparcos parallax from Ivan Leeuwe 3 J2007h . For 
HD 122563, we estimate an interstellar extinction of Ay = 0.01 
mag based on its galactic coordinates and distance. The bolo- 
metric flux Fboi is obtained by combining my, Ay, and the bolo- 
metric correction BCy, where BCy is obtained by i nterpolating 
the tables for giant stars from lAlonso et alJ (Il999bl) . We started 
with an initial r eff of 4530 K (and [Fe/H] = -2.5) to calculate 
BCy from the tables, and then used this value to determine an 
initial Fb i- Using the initial Fboi and the derived #ld we deter- 
mined T e ff (see below). This new r e fj was then used to rederive 
BCy, Fboi and T e g, and we iterated until we converged on the 
final r eff of 4582 K using BCy = -0.472 for W and 4598 K 
using BCy = -0.466 for 6^. We note that adopt i ng the se T e s 
and interpolating the tables from lHoudashelt et al.l d2000h yields 

BCy within our error ba rs (BCy 46). Fo r Gmb 1830, we 

used Fboi and Ay from iBovaiian et alJ d2012l) . and then indi- 
rectly calculated BCy. We do not subsequently use BCy in this 
work but we report the value for reference. Both the ID and 3D 
limb-darkened angular diameters #ld are given for HD 122563 
and the ID diameter is given for Gmb 1830 (see Table Rp). The 
surface brightness relations from iKerve 11a et aD (2004a c) were 
used to provide an estimate of the ID limb-darkened angular 
diameter 9 pie d. The predicted values are lower than the derived 
values, although for HD 122563 the agreement is quite good 
(#pred = 0.928 mas, 0\£> = 0.948 mas). These relations have 
been calibrated with a large sample of stars. However, the obvi- 
ous lack of reliable measurements of metal-poor stars may lead 
to slight biases in the angular diameters predicted using these 
methods. 

Combining the above mentioned measurements we deter- 
mined the observed or model- independent fundamental proper- 
ties of both stars; absolute magnitude My, T e g, L, and R, where 

/ , F \0.25 

r e ff is derived using the equation T e g — (^"gr) > °"sb is the 
Stephan-Boltzmann constant and 9 is the limb-darkened angular 
diameter. 

The most recently published NLTE spectroscopic analysis 
of HD 122563 yielded [Fe/H] = -2.53 ± 0.02 dex (see Table [TJ 



iMashonkina et al.|l2008l) . The same authors also derived NLTE 
abundances for two a elements: [Mg/H] = -2.2 and [Ca/H] = 
-2.3 to -2.4 dex. In the PASTEL catalogue there are 15 spectro- 
scopic determinations of [Fe/H] since 1990 (mostly LTE) with 
a mean value of -2.7 dex, or 5 determinations since 2000 with 
a mean of -2.6. The mean metallicity, which is a mixture of Fe 
peak and a elements then becomes [Z/X] s = -2.3 + 0.1 dex. 
Spectroscopic \ogg values typically vary between 1.1 and 1.5 
dex (see Table [TJ. 

For Gmb 1830, iGehren et all d2006l) derived [Fe/H] = 
-1.35 ± 0.10 dex from Fe II lines, which are not supposed to 
be affected by NLTE, and an NLTE [Mg/Fe] abundance of +0.3. 
Using the latter for all a elements, this implies a [Z/X] s = 
-1.3 + 0.1 dex. The spectroscopic log g of th is star has been esti- 
mated to be -4.70 dThevenin & Idiartil 19991) from NLTE studies 
but using a temperature hotter by about 200 K. The T s g reported 
in this work would result in a downward revision of this number. 



4.2. CESAM2k models 

In order to interpret the observations of HD 122563 and 
Gmb 1830 we used the CESAM2k stellar evolu tion and struc- 
ture code dMorellll997l:lMorel & Lebretonll2008h . We tested the 
models using thr ee different equations of state (EOS): the clas- 
sical EFF EOS dEggleton et ail 1 19731) with/witho ut Couloumb 
corrections (CEFF/EFF), and the OPAL EOS dRogers et all 
fl996l) . and we found small differences in the derived pa- 
rameters for Gmb 183 only. For all of the m odels we used 
the OPAL opacities dRogers & Iglesiasl 1 19921) supplemented 
with lAlexander & Ferguson! d 19941) molecular opacities. The 
p-p chain, CNO, and triple-q nuclear reac tions were calcu- 
lated using the NAC RE rates dAngulolll999l). We adopted the 
solar abundances of iGrevesse & Noelsl dl993h (Z G = 0.017, 



0.694) and u sed the MARCS model atmospheres 



4 The IRFM allows one to calculate by comparing the ratio of 
infra-red to bolometric flux observed from Ea rth to the true intr insic 
value obtained from theoretical models (see e.g.lCasagrandel l2~008h 



iGustafsson et al.l d2003l) . Microscopic diffusion was taken into 
account for Gm b 1830 and follows the treatment described by 
iBurgersid 19691) . and we introduced e xtra mixing by employing a 
parameter, Rev = 1, as prescribed by iMorel & Theveninl d2002h 
in order to slow down the depletion of helium and heavy el- 
ements. For HD 122563 no observable difference is found be- 
tween diffusion and non-diffusion models for giants, except a 
small effect on the age of the star, i.e., for the same parameters 
the model with diffusion fits the observational data with an age 
~0.3 Gyr older than the non-diffusion model. Convection in the 
outer enve lope is treated b y using the mixing-length theory de- 
scribed bv lEggletonldl972l) . where I = aH p is the mixing-length 
that tends to as the radiative/convective borders are reached, 
Hp is the pressure scale height, and a is an adjustable parameter. 
To match the solar luminosity, T e g, and oscillation frequencies 
(while including diffusion) we find a value of a — 2.04. We note 
that we did not include convective overshooting in our models 
because the primary effect that this extra parameter has on the 
determination of the stellar model is the age. This means that it 
is possible to find two equivalent stellar models with the same 
stellar parameters that differ only by age and the value of the 
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Table 7. Observed parameters of HD 122563 and Gmb 1830. 



Observation 



m v (mag) 
m K (mag) 
7T (mas) 
[Z/XL (dex) 
BCy (mag) 
Ay (mag) 

F bol (erg s^'crrr 2 x 10~ 8 ) 

^red ( mas > 
B LD (mas) 

My (mag) 
Teff (K) 
L (Lq) 
^(Ro) 



HD 122563 
ID 



Gmb 1830 

3D 

6.45 ± 0.02 
4.37 ± 0.03 
109.99 ± 0.41 
-1.3 ±0.1 

-0.23 ± o.or 

0.00 ± 0.01 
8.27 ± 0.08 
0.630 ± 0.013 
0.679 ± 0.015 



6.19 ±0.02 
3.69 ± 0.04 
4.22 ± 0.35 
-2.3 ±0.1 
-0.472 ± 0.02 
0.01 ±0.01 
13.23 ±0.37* 
0.928 ± 0.019 
0.948 ± 0.012 

-0.69 ± 0.03 
4585 ± 43 
232 ± 6 
24.1 ± 1.9 



-0.466 ± 0.02 

13.16 ±0.36* 

0.940 ±0.011'' 

-0.69 ± 0.03 
4598 ± 41 
230 ± (f 
23.9 ± 1.9 



6.66 ± 0.02 
4818 ±54 
0.213 ± 0.002 f 
0.664 ±0.015 



Notes. "BCy is derived assuming F bo i from lBovaiian et all d2012l) . b F ho \ is derived using my , Ay, and BCy from lAlonso et al.ldl999bl) . c 8 pK i is the 
predicted angular diameter using the surface-brightness relations from lKervella et afl d2004al lcl). rf 3D limb-darkened angular diameter. e L calculated 
from F bol and n. 



overshoot parameter. Since we have no observable constraint to 
distinguish between these two parameters we chose not to in- 
clude it. 

Each stellar model is defined by a set of input model pa- 
rameters — mass M, initial helium content Y\, initial metal to 
hydrogen ratio ZJX\, age t, and the mixing-length parameter a 
— and these result in model observables, such as a model T e ff 
and a model L. By varying the parameters M, Y{, Z{jX{, t, and 
a we aimed to find models that fitted the luminosity, T e s, and 
metallicity constraints as outlined in Table [7] We stopped the 
evolution of the models when an age of 14 Gyr was reached. For 
HD 122563 we chose to use the constraints from the more realis- 
tic 3D models, although we note that the difference between the 
ID and 3D constraints leads to only very slight changes in the 
parameters of the stellar models (see Sect. l4.3.Tl below). 

4.3. Stellar parameters 

Figure [4] shows two HR diagrams with the observational error 
boxes of both stars (left/right = HD 122563/Gmb 1830) as well 
some models that lie somewhat away from the central position of 
the box, illustrative of the uncertainties that we find in the stellar 
parameters (see below). Table [8] lists the stellar parameters for 
both stars using the classical EFF, and for Gmb 1830 we also 
give the stellar properties for the CEFF and OPAL EOS models. 

Given the few independent observational constraints and the 
large number of adjustable parameters in the models, a classical 
error analysis is not possible for both stars. In order to estimate 
the uncorrelated uncertainties we changed each of the reference 
parameters of the models individually until we reached the edges 
of the error box in the HR diagram, or the limits of each param- 
eter, e.g. we did not test Y\ < 0.20. These are the uncertainties 
that are given in the top part of Table [8] For the uncertainty in 
the age we give the 1 cr uncertainty which corresponds to the cen- 
tral models approaching the upper and lower limit in luminosity 
(first number) and we also give the range of possible ages while 
considering the uncertainties in the four model parameters (sec- 
ond uncertainty). We also list the model observables and their 
uncertainties in the lower part of the table. We note that the un- 
certainties in the model observables cover the full range of val- 



ues while considering the individual changes in each of the four 
model parameters. 

4.3.1. HD 122563 

For HD 122563 using the EFF EOS description we found a best 
model with M = 0.855 M , Y, = 0.245, a = 1.31, and t = 12.6 
Gyr. We fixed Z-JXi = 0.0001 in order to have the correct ob- 
served [Z/X] s . This model is illustrated in Figure |4] (left panel) 
by the thick line and is clearly labelled. We also show some mod- 
els which illustrate the reported uncertainties: the black contin- 
uous line shows a model when a is changed by 0.08, the red 
continuous line shows the central model when Yi is decreased 
by lcr, and finally the red dashed-dotted line shows the effect 
of increasing the mass by lcr. We note that if we increase the 
mass to more than 0.88 M then the age of the model becomes 
too small (< 10 Gyr) if we are to consider the giant a halo star. 
We also found correlations among the parameters M, Y{, and a, 
and adjusting two of the three at a time by a small amount re- 
produces the position of the central model, e.g. if we fix M then 
AY = +0.01 <=> Aa = -0.01. However, these correlations are 
adequately accounted for in the uncertainties. 

The dotted error box shows the constraints if we consider 
the ID limb-darkened angular diameter. The stellar parameters 
of the model that passes through the center of the box need small 
adjustments when compared to the 3D diameter constraints. In 
particular, decreasing either M or Y, alone by less than lcr or 
decreasing a by ~ 0.03 (or a combination of the three) would 
reproduce the central position of the error box with a slightly 
higher age. If the temperature constraint were even lower, then 
the only viable option would be reducing the mixing-length pa- 
rameter a, because reducing M or F,- by much more would result 
in a model that fails to reach the minimum luminosity before 14 
Gyr. 

Inspecting the stellar parameters in Table [8] we highlight 
the excellent precision obtained in the mass of this single star. 
Generally such precisions can only be obtained if the star is in a 
binary system, where the solutions are then model-independent. 
Combining this value with the well-determined radius yields a 
very precisely determined logg (= 1.60 + 0.04 dex). This value 
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is larger than most values used for spectroscopic analyses which 
typically ranges from 1.1 - 1.5 dex (see Table [TJ. More recent 
work using 3D hydro-dynamical simulatio ns for stellar atmo- 
spheres quote values of 1 . 1 - 1.6 dex (se e e.g. lBarbuv et al.ll2003l 
ICollet et al. 2001 iRamfrez et al.ll201O) . 



4.3.2. Gmb1830 

In Table |8]right three columns we summarize the stellar param- 
eters for Gmb 1830 using the EFF, CEFF, and OPAL EOS. In 
Figure [4] we show the central model for the EFF EOS (arrow 
with 'EFF') with illustrative uncertainties. The model parame- 
ters are M = 0.635 + 0.025 M , Y, = 0.235 ± 0.025, [Z/X] = 
0.0016 ± 0.0004, a = 0.68 + 0.10, and t = 12.0 + 02 + _\\. We 
also show a CEFF and OPAL EOS evolution track using the cen- 
tral parameters obtained with the EFF model. A qualitative dif- 
ference between the three EOS is notable, however, considering 
the uncertainties in the stellar parameters, these differences are 
not significant. 

The uncertainties reported in Table [8] do not consider all of 
the correlations among the parameters. For example, reducing 
the mass by lcr implies a necessary increase in T; by lcr in or- 
der to remain inside the error box and vice versa. In Figure |4] we 
show effects of the uncertainties on the central model; the dotted 
black line shows the effect of decreasing a by 0.10, the dashed 
black line shows the effect of decreasing Z\jX\ by lcr (denoted 
by AZ in figure), and the red continuous line right of the central 
model is when the mass is decreased by lcr and Y\ increased by 
lcr. We note that by decreasing/increasing the mass or Y\ alone 
leads to a very young stellar model (not consistent with a halo 
star), a model that is too hot, or at the age of 14 Gyr the lumi- 
nosity does not reach the minimum required 0.210 L . 



4.4. Asteroseismic Constraints 

In Table [8] we predict two global ast eroseismic quantities (Ay) 
and Vm ax based on scaling relations (iBrown & Gillilandl[l994l 
iKieldsen & Beddingll 19951) corresponding to the reference mod- 
els. Both quantities are proportional to the mass and radius of 
the star, with the latter also having a small r e ff-dependence; 



<AV> ~ 



M 



<Av> G 



R 2 Vr eff /5777K 



(1) 



where (Ay) 



134.9 uliz and 



„„ t>0 = 3,050 //Hz 
dKjeldsen & Beddingll 19951) . and R and M are in solar units. 
Although the relations are approximate, they have been found 
to wor k quite well, e.g. lBedding & KieldsenM2003l) : IStello et alj 
(120081) . (Av) is the characteristic spacing between consecutive 
radial overtones of the same mode degree seen in the power 
(frequency) spe ctrum of a star w ith sun-like oscillations (e.g. 
see Fig. 6 from lButler et al.ll2004l) . and it is proportional to the 
square root of the mean density of the star. Because it is a repeti- 
tive pattern (similar to a periodicity), it is relat ively easy to deter- 
mine from even low signal/nois e data (see e .g. Huber et al.l2009t 
Mosser & Appourchauxl 120091; iRoxburghl 120091; iMafhur et al.l 
2010t IVerner et al.l 1201 ll who discuss different methods to de- 
termine this value). The value of v max is the frequency corre- 
sponding to the maximum amplitude of the bell-shaped ampli- 
tude/power spectrum, and it is also a quantity that can be more 
easily observed than, for example, individual oscillation modes. 

Because the radii and effective temperatures of these stars are 
well determined, the predicted seismic quantities depend only 
on the mass of the star. If we substitute directly the derived 



mass ranges into the equations then we can predict the range of 
possible values for these quantities corresponding to the central 
model, i.e. not taking into account the changes in a, Y{ or Z\jX\. 
For Gmb 1830 we find that for masses = [0.62, 0.64, 0.66] M Q 
we calculate v max = [4773, 4927, 5081] /iHz and <Av) = [196, 
199, 203] /iHz, which correspond to typical periods of approxi- 
mately 4 minutes. If we can detect these values, even with poor 
precision we will still be able to select the optimal mass range 
and discard certain solutions. We note that both M and Y\ are 
very highly correlated, and so fixing M will present interesting 
constraints on Y\. Performing such observations from ground- 
based instrumentation should yield successful results (e.g. a 2+ 
meter telescope equipped with a highly efficient and stable spec- 
trograph). ForHD 122563 we find forM = [0.84, 0.86, 0.88] M Q , 
v max = [5.03, 5.15, 5.26] pHz and (Av) = [1.05, 1.06, 1.07] ^uHz. 
The dominant periods are approximately 2.5 days, and observa- 
tions from ground-based instrumentation would be difficult. In 
order to use asteroseismic data to help constrain the models for 
HD 122563, we would require seismic data from space-borne 
instruments, such as with the CoRoT or Kepler missions, to pro- 
vide the necessary precision and determine the individual oscil- 
lation modes. 



5. Conclusions 

We have determined the T e g, L, and R of HD 122563 and 
Gmb 1830 by using K band interferometric measurements 
(Table [7} and 3D/1D limb-darkening for the giant/dwarf. We 
find angular diameters of 0t, d = 0.940 ±0.011 mas and 9\d = 
0.679 + 0.015 mas for HD 122563 and Gmb 1830, respectively, 
and these convert into T eff = 4598 + 41 K for HD 122563 and 
T e g = 48 1 8 + 54 K for Gmb 1830. These new precision tempera- 
tures increase the well-known difficulty of fitting the error boxes 
of these two metal-poor stars with evolutionary tracks. Using the 
CESAM2k stellar structure and evolution code we found that we 
could match models to the data by using values of the mixing 
length (the parameter a) very different from that of the Sun. We 
found values of a = 0.68 and 1.31 for the 0.63 M dwarf star 
and the 0.86 M giant, respectively. The or der of these values 
seems consistent wit h recent model analyses dYildiz et al.ll2.006h 
iKervella et alJl200~8h . We found that different equations of state 
lead to qualitatively but not quantitively different model param- 
eters for the dwarf star but not for the giant. The initial helium 
content comes out similar to the big-bang value, the deduced 
masses are low and their ages are high, consistent with expected 
values for metal-poor halo stars (see Table [8]). The masses are 
determined with a few percent precision and coupling these with 
the radii yields well-constrained values of log g. For the giant 
star we found logg = 1.60 + 0.04 somewhat higher than the 
typical values (1.1 - 1.5) adopted by spectroscopic anal yses ac- 
cording to the PASTEL catalogue ( ISoubiran et al.ll2010|) and for 
the dw arf star we obtain \ogg = 4.59 ± 0.02 dex. iBarbuv et alJ 
(120031) determined the O abundance of HD 122563 assuming 
two different (both justified) values of log g, and they concluded 
that their resulting [O/Fe] = +0.7 abundance seemed most con- 
sistent when they adopt the Hipparcoffl log g — 1 .5 and not the 
value determined from ionization equilibrium of Fe, logg = 
1.1, a result due possibly to NLTE effects. This work supports 
their O determination. With both \ogg and r e ff now very pre- 
cisely known, these provide very important inputs for any spec- 
troscopic analyses, especially for the determination of neutron- 



5 We note that with the new Hipparcos parallaxes the deduced 
logg£ 1.6. 
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Fig. 4. HR diagram showing the observational error boxes for HD 122563 (left) and Gmb 1830 (right). Both figures show stellar 
models that pass through the error boxes which allow us to determine the stellar model properties and their uncertainties. Each 
panel shows the adopted central models (with arrows) obtained by considering the HR and metallicity constraints. Other models are 
also indicated to highlight the parameter uncertainties and correlations. Refer to Sects. 14370 and l4~3.2l for details. 



Table 8. Stellar properties and lcr uncertainties derived from modelling HD 122563 (no atomic diffusion) and Gmb 1830 (with 
atomic diffusion). 





HD 122563 Gmb 1830 






Property 


EFF EOS 


CEFF EOS 


OPAL EOS 



M(M ) 0.855 + 0.025 0.635 ± 0.025 0.625 ± 0.015 0.620 + 0.020 

Y, 0.245 + 0.015 0.235 + 0.025 0.230 + 0.020 0.235 + 0.025 

Zj/X, 0.00010 + 0.00002 0.0016 + 0.0004 0.0016 + 0.0004 0.0016 + 0.0004 

a 1.31 +0.08 0.68 + 0.10 0.63 + 0.08 0.65 + 0.10 

Age(Gyr) 12.6±0.1!|° 12.1 + 0.2+^ 12.7 + 0.3^ 12.3 ± 0.3^ j 

R(R e ) 24.1 + 1.1 0.665 + 0.014 0.665 + 0.015 0.665 + 0.015 

L (L e ) 230+ 7 0.213 + 0.002 0.213 + 0.002 0.213 + 0.002 

T eff (K) 4598 + 42 4815 + 52 4814 + 53 4815 + 50 

logg(dex) 1.60 + 0.04 4.60 + 0.02 4.59 + 0.02 4.58 + 0.02 

[Z/X] s -2.38 + 0.10 -1.32 + 0.11 -1.32 + 0.11 -1.33 + 0.11 

<Av>° rcd OuHz) 1.06 + 0.06 198 + 6 197 + 7 196 + 6 

5.16 ± 0.38 4886 ± 190 4809 ± 199 4768 ± 188 



Notes: The first five values are the input parameters of the model and the other values are properties of these models. The uncertainties are derived 

by perturbing each of the model parameters individually until the edge of the error box is reached. 

°{Av) and v max are the predicted seismic quantities according to the scaling relations from lBrown & GillilanJ ( fT99 4): Kieldsen & Bedd inddl995l) . 
and the range of values listed consider all the uncertainties in the five model parameters. 



capture element abundances which can constrain models of nu- 
cleosynthesis. 

Finally, we have also predicted the asteroseismic signatures 
(Av) and v max for the two stars and we showed that determi- 
nations of these quantities for the dwarf star are possible using 
ground-based observations. For the giant, however, we would 
require very long time series in order to resolve the frequency 
content of the oscillations, and this would only be possible with 
space-borne instruments. The asteroseismic data would provide 
very important constraints because it would allow us to deter- 
mine the mass with better precision (using the radius from this 
work), and thus the initial helium abundance. 
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